Spatially improved operators for excited hadrons on the lattice 
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We present a new approach for determining spatially optimized operators that can be used for 
lattice spectroscopy of excited hadrons. Jacobi smeared quark sources with different widths are 
combined to construct hadron operators with different spatial wave functions. We use the variational 
method to determine those linear combinations of operators that have optimal overlap with ground 
and excited states. The details of the new approach are discussed and we demonstrate the power of 
the method using examples from quenched baryon and meson spectroscopy. In particular we study 
the Roper state and p(1450) and discuss some physical implications of our tests. 
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I. INTRODUCTION 

Ground state spectroscopy on the lattice is by now a 
well understood physical problem and impressive agree- 
ment with experiment has been achieved. The lattice 
study of excited states is not so far advanced and for 
many systems even the correct level ordering of the states 
has not yet been observed. A prominent example for such 
a system is the nucleon and its excited states. In particu- 
lar the first excited positive parity nucleon, the so-called 
Roper state, has seen a lot of attention from the lattice 
community during the last two years 0-3- Revealing 
the true nature of the Roper state with non-perturbative 
methods is an important task. 

In a lattice calculation the masses of excited states 
show up in the subleading exponentials of Euclidean two 
point functions. A direct fit of a single Euclidean correla- 
tor is cumbersome since the signal is strongly dominated 
by the ground state. Also with methods such as con- 
strained fits [j| or the maximum entropy method [|| one 
still needs very high statistics for reliable results. 

An alternative method is the computation of not only 
the correlator of a single operator, but the calculation 
of a full matrix containing all cross correlations of a set 
of several operators with the correct quantum numbers 
0- In this so-called variational method the correlation 
matrix is then diagonalized and one can show that each 
eigenmode is dominated by a different physical state (for 
a properly chosen set of basis operators). After normal- 
ization at distance 0, for t > the largest eigenvalue 
gives the Euclidean correlator of the ground state, the 
second-largest eigenvalue corresponds to the first excited 
state, et cetera. 

The success of the variational method depends strongly 
on the choice of the basis operators. They need to be lin- 
early independent and should give rise to a large overlap 
with the physical states. A physical hadron state has sev- 
eral characteristics, in particular a specific Dirac struc- 
ture and its spatial wave function. Therefore it is advan- 
tageous to optimize the spatial properties of the interpo- 
lating operators. An example for this fact is the men- 



tioned Roper state where the variational method based 
on nucleon operators that differ only in their diquark con- 
tent but have the same spatial wave function did not lead 
to success ,4] . It can be argued that a node in the radial 
wave function is necessary to capture reliably the Roper 
state or other radially excited hadrons. 

In lattice calculations it is important to devise an effi- 
cient way of implementing a spatial wave function at low 
numerical cost. When the quark wave function extends 
over several lattice sites, a naive approach would require 
the inversion of the Dirac operator on point sources lo- 
cated at all possible quark positions (see 8] for a discus- 
sion of such wave functions) . Creating a wave function for 
quarks in this way quickly becomes numerically expen- 
sive. For improving the ground state wave function an 
effective technique, so-called Jacobi smearing, has been 
shown to give good results @- Here, a point-like quark 
source is smeared to a shape similar to a Gaussian, i.e. 
the origin of the source is connected to neighboring lattice 
sites within a time slice by gauge transporters. Such a 
source increases the overlap of the physical hadron states 
with the lattice operators used to create these states and 
considerably reduces the fluctuations of hadron propaga- 
tors. 

In this article we demonstrate that combining Jacobi 
smeared quark sources with different widths in the varia- 
tional method provides a powerful tool for the analysis of 
excited hadron states. After presenting the outlined ideas 
in detail we apply our method to quenched baryon and 
meson spectroscopy. The excited nucleon system as well 
as excited mesons are analyzed (for recent lattice stud- 
ies of excited mesons see [l(| )• We find good effective 
mass plateaus for the first and partly the second radially 
excited states. The propagators can then be fitted us- 
ing standard techniques. Different physical ramifications 
and implications of our findings are briefly discussed. 

II. THE METHOD 

Let us begin the presentation of our method with a 
brief recapitulation of Jacobi smearing of quark sources. 
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Since a complete quark propagator (the inverse of the 
lattice Dirac operator D) is a far too large object to 
be stored completely in the computer, one has to work 
with the propagator evaluated on some source s^ a ' c \ The 
source is placed at timcslicc t = of the lattice. It is la- 
beled by a Dirac index a and a color index a. One then 
computes {p and c are summed over) 



Narrow source 



Wide source 



„(a>°) 



(x,t) = ^D^&t | y,0), P s (Q ' a) (y,0)p . 



(1) 

When one chooses a point-like source at the spatial origin 
0, i.e. s — sq with 



4 Q,< %0)p = 6(y~0)S pa 6 ca 



(2) 



the resulting vector is the quark propagator from the 
origin to all lattice points. These quark propagators can 
then be combined to form hadron propagators. 

Choosing a point-like quark source has, however, the 
big disadvantage of a poor overlap with the true wave 
function. After all, for the lattice spacings which are 
appropriate, quarks are not expected to be located at 
a single point inside the hadron and the overlap of the 
point-like wave function with the true physical wave func- 
tion is small. The situation can be improved, e.g. by Ja- 
cobi smearing. One acts with a smearing operator M on 
the point-like source Sq to obtain the smeared source 

s (a,a). 
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Jf J Q > a ) 
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H(x,y) = (a-, 0)«5( 



+ Uj(x-j ,0)U{x-j ,y) 



(3) 



The operator H is simply the spatial hopping part of the 
Wilson term at timeslice 0. Note that H, and thus M, is 
trivial in Dirac space and acts only on the color indices 
(in our notation we suppress the color indices of M, H 
and of the gauge transporters U). 

The Jacobi smearing outlined in Eq. © has two free 
parameters: The number of smearing steps N and the 
positive real parameter k. These two parameters can be 
used to adjust the profile of the source. In order to study 
the profile of the smeared source s^ '^ we define 



P(r) = 



y 



E 



= ( Q . a ) 



(y,o)< 



(4) 



Since the smearing is trivial in Dirac space the profile 
function P(r) is independent of the Dirac label a which 
appears in the source s^ a ' a \ The profile function P(r) 
can, however, depend on the color label a, but as we will 
demonstrate below this dependence is small. The delta 
function on the right hand side is implemented by binning 




FIG. 1: Profiles P(r) of the narrow and wide source. The 
symbols are our data points, the curves are the target Gaus- 
sian distributions which we approximate by the profiles P(r). 



the discrete values of \y\ on the lattice. We stress that 
P(r) is not a gauge invariant object, but only an auxiliary 
quantity to visualize the source. Certainly it is easy to 
construct a profile function which is a color singlet, but 
since we are interested in visualizing the sources s( Q > a ' 
individually for all a = 1, 2, 3, we use the form Eq. (JIJ. 

In this article we work with two different sources, a 
narrow source n and a wide source w with parameters 
given by 



narrow source n 
wide source w 



N = 18, k = 0.210 

N = 41 , k = 0.191 (5) 



In Fig. we show the profile functions P(r) for these 
two sources. They were normalized such that P(0) = 1. 
For two different gauge configurations (20 3 x 32, Liischer- 
Weisz gauge action, (3 = 7.90, lattice spacing a = 0.148 
fm) we superimpose the values for P(r) for all three pos- 
sible color indices of the original point source, i.e. for 
a = 1,2,3. It is obvious from the plots that the differ- 
ent color components of the point source lead to very 
similar profiles for the smeared source. Also when com- 
paring different configurations, the relative fluctuations 
are small. The smearing parameters N and k were cho- 
sen such that the profiles approximate Gaussian distri- 
butions with widths of a ~ 0.27 fm for the narrow source 
and a ~ 0.41 fm for the wide source. These Gaussians 
are displayed as curves in Fig. ^ 

We remark that the parameters were chosen such that 
simple linear combinations c n n + c w w of the narrow 
and wide profile approximate the first and second ra- 
dial wave functions of the spherical harmonic oscillator: 
The coefficients c„ ~ 0.6, c w ~ 0.4 approximate a Gaus- 
sian with a width of a ~ 0.33 fm, while the combination 
c n ~ 2.2, c w ~ —1.2 approximates the corresponding ex- 
cited radial wave function with one node. Thus, the two 
sources we include allow the system to build up radial 
wave functions with and without a node. 
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The final form of the wave function is, however, not 
put in by hand, but is determined by the system through 
the variational method Q ■ In this approach one does not 
calculate a single correlator, but instead a complete cor- 
relation matrix of operators Oi,i — 1,2, ... R that create 
from the vacuum the state which one wants to analyze. 
One calculates all cross correlations 

C(% = (O,(t)O}(0)) . (6) 

In Hilbert space this correlation matrix has the represen- 
tation (for infinite temporal extent) 

C(% = ^(0|O l |^)(«|O]|0) e -* M " , (7) 
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where the sum runs over all physical states |n) and the 
corresponding energies are denoted as M n . The eigen- 
values A^" (t) of the correlation matrix can be shown to 
behave as 

A (fe) (i) oc e -t Mk [ 1 + 0(e -t AMfc ) ] , (8) 

where AM& is the distance of Mk to nearby energy lev- 
els. A modification of the method is the analysis of the 
generalized eigenvalue problem 

C(t)v = \{t)C(t )v, (9) 

which can be rewritten to a standard eigenvalue problem 
by bringing C(to) to the left hand side of Eq. ©, e.g. in 
the form of C(io) -1 multiplied from the left. The nor- 
malization at some slice to < t is expected to improve the 
signal by suppressing the contributions of higher excited 
states. We will explore the freedom of such a normaliza- 
tion in our analysis. 

The sources we use for the correlation matrix are con- 
structed from the narrow and wide quark sources we pre- 
pared. This is best explained in an example: An operator 
which creates the p-meson is given by ujid. Both the u 
and the d quark can either have a narrow (n) or a wide 
(w) quark source. This gives the four possible combi- 
nations (n, n), (n, to), (to, n), (w, w), where the first entry 
refers to the smearing of the u quark and the second en- 
try is for the d quark. Thus, we can use a basis of the 
four operators 

Oi = (n,n), 2 = (n,w), 3 = (w,n), C 4 = (to, to) , 

( 10 ) 

for building up the correlation matrix C(i). We remark 
that the different operators Oj also have to be used at 
the sink end (timeslice t). This can be implemented by 
applying the smearing operator M with the two sets of 
smearing parameters Eq. JSJ at the sink end of the quark 
propagator. 

Before we put our method to a test in the nucleon 
and meson systems let us briefly summarize some tech- 
nical details. For our quenched calculation we use the 
chirally improved Dirac operator It is an approxi- 
mation of a solution of the Ginsparg Wilson equation ^1 



which governs chiral symmetry on the lattice. The chi- 
rally improved Dirac operator is well tested in quenched 
ground state spectroscopy where pion masses down 
to 250 MeV can be reached at a considerably smaller 
numerical cost than needed for exact Ginsparg Wilson 
fermions. For ground states the chirally improved action 
shows very good scaling behavior. 

The gauge configurations we use for testing our method 
were generated on a 12 3 x 24 lattice with the Liischer- 
Weisz action The inverse gauge coupling is f3 = 7.9, 
giving rise to a lattice spacing of a = 0.148(2) fm as 
determined from the Sommer parameter in 15]. The 
statistics of our ensemble is 100 configurations. We use 
10 different quark masses m ranging from am = 0.02 to 
am = 0.20. 



III. EXAMPLE A: EXCITED NUCLEONS 

The first example where we put our new method to a 
test is the spectroscopy of excited nucleons. Spectroscopy 
of the lowest positive and negative parity states is a key 
to understanding the physics of the nucleon. From a 
lattice perspective the nucleon system still holds a few 
unresolved puzzles and new methods will help to obtain 
a more complete picture. 

It has long been noted that the observed ordering of 
the lowest positive, l/2 + , iV(1440), and negative parity 
excitations of the nucleon, l/2~, iV(1535) is 'unnatural'. 
Indeed, a physical picture based on linear confinement, 
Coulomb and color-magnetic terms, always arranges the 
first radial excitation above the first orbital excitation, 
i.e. the excited states have alternating parities. This 
outcome is in contrast to the observed nucleon masses, 
where the first radial excitation, AT(1440), is well below 
the first orbital one, iV(1535). This has prompted spec- 
ulations that perhaps the Roper resonance is not a three 
quark state, but a collective excitation of the bag surface 
m|, a gluonic state 01 > a No coupled channel effect 
[l8j , a resonance in the pion - skyrmion system , and 
most recently - a pentaquark state with a scalar diquark 
- scalar diquark - antiquark structure |20j . 

On the other hand, we have to expect that close to the 
chiral limit effects of the spontaneous breaking of chiral 
symmetry should be important. It has been suggested 
in |2l| that in the low-lying baryons the residual interac- 
tion between the valence constituent quarks mediated by 
the Goldstone boson field should be of vital importance. 
Such an interaction is of the flavor- and spin-exchange 
nature, contrary to the perturbative QCD degrees of free- 
dom, and very naturally resolves the puzzle of low-lying 
baryon spectroscopy in the u,d,s sector. 

If chiral symmetry breaking is important for the Roper 
state, we expect that the ordering of the lowest excita- 
tions of baryons should change as a function of the quark 
mass. Baryons made of heavy quarks, where chiral effects 
do not play a role, should show a radial excitation above 
the orbital excitation, while for smaller quark masses we 
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expect a level reordering as seen in the nucleon system. 
At intermediate quark masses a level crossing of the ra- 
dial and orbital excitations should take place. Hence, 
studying the evolution of the baryon spectrum versus the 
current quark mass allows one to clarify the physical pic- 
ture. This can be done within the lattice approach where 
masses of quarks are external parameters which can be 
freely varied. In particular with the newly developed 
fermion actions based on the Ginsparg- Wilson equation 
the region of small quark masses also becomes accessible. 

Previous attempts to study these issues on the lattice 
Q - 3 have mainly used the two interpolating fields, xi 
and X2, 

Xl(x) = tabc [Ua( X ) C75 d b {x)] U c (x) , (11) 

X2(x) = e abc [v£(x)C d b (x)] 7s u c (x) . (12) 

Here C is the charge conjugation matrix and a, 6, c are 
the color indices. For the Dirac indices we use matrix no- 
tation. While the coupling of the operator xi to both the 
nucleon and the lowest negative parity state has been re- 
liably established, the operator X2 couples neither to the 
nucleon, nor to the Roper state. For a detailed study of 
this issue and for an interpretation of this fact see 0. 
Hence a remaining possibility to sec the Roper state is to 
use Xi an d try to separate the strong signal of the ground 
state (the nucleon) from the weak signal of the radial ex- 
cited state (the Roper state, if indeed the Roper state is 
a 3-quark state) at small Euclidean time. This strategy 
has been followed in Refs. and where in the former 
case a multi-exponential fit of the diagonal (xiXi) corre- 
lator has been performed (using constrained curve fitting 
la), while in the latter the maximum entropy method 
01 has been used. Both papers report the observation of 
the Roper resonance and level crossing towards the chiral 
limit. However, the two papers contradict each other on 
the value of the pion mass where this level switching takes 
place (300-400 MeV versus 600 Mev). Further improve- 
ment of the lattice technology is necessary to deepen the 
understanding of the nucleon system from the lattice. 

Both the open physical and technical questions of the 
nucleon system make it an ideal testing ground for our 
new approach. The goal is to see whether this method 
allows one to identify a reliable signal for the excited 
positive parity nucleon. Our analysis is based on the 
interpolator Xi defined in l|ll|) . It contains three quarks 
and each of these quarks can be smeared either narrow 
(n) or wide (w). This gives 8 possible combinations 

0\ = (n,n,n), O2 — (n,n,w), ... Os — (w,w,w) . (13) 

In this notation the first entry refers to the left-most 
quark in Eq. (|11|) . i.e. the u quark inside the diquark 
part of xi- We remark that 'diquark' does not refer to 
a true clustering of quarks in the nucleon but is used 
for the combination of the first two quark operators in 
the interpolators Xi?X2- The second entry is for the d 
quark and the third entry for the other u quark outside 



the diquark in xi- From these operators we calculate the 
correlation matrix 

Cg(t) = (Oi(t) ~[1 ±74]Oi(0)> , (14) 

where we have inserted projectors to positive and nega- 
tive parities. Using the relation C + (t) — — C~(T — t), 
where T is the total time extent of our lattice, we com- 
bine the two correlators to improve the statistics. This 
gives rise to the combined correlator C(t) which we then 
use in the variational method. The signal for positive 
parity states is obtained for small t when running for- 
ward in time, while the negative parity states propagate 
backward in time (T — t). Near T/2 there is a cross- 
ing region where the propagators of the two parities mix. 
We focus on the positive parity states and show plots of 
propagators and effective masses only up to about T/2. 

The correlation matrices C{t) are real and symmetric 
within error bars and we symmetrize the matrices by re- 
placing Cij(t) with [dj (f) + Cji(t)]/2. Subsequently, we 
calculate the eigenvalues \( k '(t) for all t. At each time- 
slice t we order the (real) eigenvalues, such that A^(i) is 
the largest eigenvalue, A^ 2 ^ (t) is the second largest eigen- 
value et cetera. We remark that it has been suggested , 
that one should analyze the eigenvalues of the generalized 
eigenvalue problem Eq. ©, e.g. by studying the eigen- 
values of the normalized correlation matrix C(to)~ 1 C(t) 
for some to < t. This normalization is expected to reduce 
the uncertainties due to the admixture of higher excita- 
tions. We do not confirm such an improvement of the 
signal. We will discuss this normalization below and for 
now continue the discussion using eigenvalues from the 
unnormalized matrix C(t). 

An important issue is the choice of the operators that 
are included in the analysis. Too few operators may not 
be sufficient to span all physical states we want to ana- 
lyze. Too many operators drive up the numerical cost 
without necessarily improving the signal. If a newly 
added operator creates a state which only has small over- 
lap with true physical states then it will not contribute to 
improving the signal. On the contrary, it will contribute 
noise to the correlation matrix and even decrease the 
quality of the signal. When exploring our new method 
we analyzed correlation matrices starting from size 2x2 
all the way up to the maximal size of 8 x 8. We find that 
the results for the masses of the ground state and the first 
excited state agree within one standard deviation (s.d.) 
when comparing different combinations of operators. We 
observe that when using more than 4 basis operators the 
signal for the lowest two states does not improve any 
further. 

For correlation matrices of four operators we systemat- 
ically analyzed the possible combinations. For some sets 
of operators, such as e.g. (n, n, n), (n,n,w), (n,w,n), 
(w, n, n) we found that the effective masses from the first 
and second excited states are degenerate within error 
bars. Indeed, the operators (n,w,n) and (w,n,n) must 
equally couple to the same physical state. This is because 
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the two quarks within the brackets of form a scalar- 
isoscalar diquark and can always be interchanged. Hence 
an unnecessary repetition of equivalent operators, like 
(n, w, n) and (w, n, n) should be avoided. When working 
with the set 

(n,n, n), (n,n,w), (n,w,n), (w, w, n) , (15) 

which does not contain equivalent operators, we find, at 
least for large quark masses, a splitting of the effective 
masses from the second and third eigenvalues larger than 
one s.d. 

To summarize the comparison of different sets of oper- 
ators we find that the masses of the ground and excited 
states agree within error bars. A distinction of the masses 
from the second and third eigenvalues is possible only for 
particular sets of operators which typically contain also 
operators with two wide sources (e.g. the set in Eq. i|15fl ). 

In Fig. we show the positive parity parts of the three 
largest eigenvalues \^(t), \( 2 \t) and X^(t) from the 
4x4 correlation matrix with the set of operators l|15|) con- 
structed from the local xi interpolator. The exponential 
decay of all three eigenvalues is clearly seen and the slopes 
differ, most obviously for the ground state and the first 
excited state. We identify these signals with the nucleon, 
the Roper state and the next positive parity resonance 
iV(1710). We remark that in the heavy quark region the 
latter two states belong to the same shell and hence must 
be approximately degenerate. However, given the fact 
that the masses of the excited states extracted from the 
second and the third eigenvalues are rather close, we per- 
form an additional test that these eigenvalues represent 
different states. 

It has been understood long ago that all Roper states 
form an excited 56-plet of SU(6). In this multiplet 
the parity of any two-quark subsystem is positive. The 
iV(1710) belongs to another multiplet, which contains 
both positive and negative parity two-quark subsystems. 
The two-quark subsystem in the brackets in xii Eq. (| X X p . 
has positive parity, while the subsystem in the brackets 
of X2 has negative parity. Hence, applying the same type 
of smearings to the xi interpolator, we should see only 
the signal from the iV(1710) state, and no signals from 
the nucleon and the Roper. We have diagonalized a 4 x 4 
correlation matrix with the source combinations listed 
in (|15l) applied to the local xi interpolator and indeed 
observed a signal only in one eigenvalue. This signal is 
clearly compatible with the 7V(1710) signal obtained with 
Xi- 

Before we discuss the effective masses and the fit pro- 
cedures we used, let us address another important advan- 
tage of our source technique. In Q evidence for nucleon- 
t]' ghost contributions (a quenching artifact) to the nu- 
cleon correlators at small quark masses was given. These 
contributions come with a negative coefficient and in 
had to be included explicitly in the multi-exponential fit- 
ting function of the correlator. In our correlation matrix 
analysis we find that the first and second eigenvalues are 
not at all affected by the ghost contributions, i.e. these 
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FIG. 2: The three largest eigenvalues A' 1 ' (circles), A^ 2 ' 
(squares) and A^ (triangles) for the quark masses am — 0.05 
(top), am = 0.10 (middle) and am — 0.20 (bottom). 



correlators are positive and we observe undisturbed effec- 
tive mass plateaus. Only the third and higher eigenvalues 
can become negative for large t. Thus, the ghost contri- 
butions are disentangled from the ground state and the 
first excited state, i.e. the states we are interested in here. 

Let us now present the effective masses as obtained 
from the correlation of the operators listed in Eq. 115fl . 
In Fig. El we show effective masses 



1 



m eff { t+ 2 ] = 111 



A(*)(t + i; 



(16) 



for the first three eigenvalues (out of a total of four) with 
the data for A' 1 ) in the left hand side column, A( 2 ) in the 
center column and A^ 3 ** on the right hand side. We display 
the data for different quark masses from am = 0.05 up 
to am = 0.20 (top to bottom). The symbols show our 
numbers for the effective masses with statistical errors 
determined with the jackknifc method. The horizontal 
lines in the plots represent our fit results: We plot the 
one standard deviation (s.d.) error band for the fitted 
mass. The band extends over the i-interval which was 
used for the fit. We stress that these are not fits to the 
effective mass, but fully correlated fits to the propagators 
shown in Fig. We discuss the details below. 

The signal for the ground state appears in the effective 
mass for A^ 1 -* which we show in the left hand side column. 
We find well pronounced plateaus. The fit range was 
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FIG. 3: Effective masses from the first three eigenvalues (left to right) for quark masses between am = 0.05 and am = 0.20 (top 
to bottom). The symbols represent the numerical data. The horizontal lines extend over the chosen fit range and represent the 
fit results plus and minus the statistical error. 



chosen from t = 4 to t = 9. For t < 4 we observe 
contributions from higher excited states, while for t > 9 
one enters the crossing region where contributions from 
the negative parity states, which propagate backwards in 
time, start to mix with the positive parity signal. 

The first excited state shows up in (center col- 
umn). It is obvious, that here the statistical errors are 
larger than for the ground state, in particular for t > 4. 
This, however, is as expected for a heavier state which 
has a faster decreasing correlator, giving rise to a smaller 
signal to noise ratio already at not too large t. However, 
also for the first excited state we find good plateaus in 
the effective mass. These plateaus extend from t = 3 
to t = 9 for our largest mass am = 0.20 and shrink to 
t = 4 to t = 7 for am = 0.05. We emphasize that there 
is a common interval in t where all 3 plateaus are seen 



simultaneously. It is very satisfactory to see a credible 
effective mass plateau for the excited states and not to 
have to rely on the fits alone to extract their masses. 

When inspecting the effective masses from the third 
eigenvalue one finds again good effective mass plateaus. 
The statistical errors are not larger than for the second 
eigenvalue and fits in similar i-ranges can be performed. 
We did not fit the A^ 3 ^ data for quark masses below 
am = 0.08 where we no longer can exclude contributions 
from ghost states and observe a decreasing quality in the 
effective mass plateaus. When comparing the positions 
for the plateaus in the first and second excited states we 
find that they are rather close to each other. Only for 
the three largest masses, where the statistical errors are 
smaller, can we establish a mass splitting larger than one 
s.d. Also in nature the mass splitting between the first 
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FIG. 4: Baryon masses as a function of the quark mass. The 
open symbols represent our lattice results and the filled sym- 
bols show the experimental masses of the corresponding nu- 
cleons, converted into lattice units with a — 0.148 fm as de- 
termined from the Sommer parameter. 



and second excited positive parity states (iV(1440) and 
iV(1710)) is relatively small. 

Before we come to presenting the final mass values of 
our analysis we need to discuss the procedures we used 
for fitting the correlators. The exponential decay of the 
eigenvalues was fitted with the two parameter ansatz 
A exp( — Mt ). Since the different values of t are not sta- 
tistically independent we used fully correlated fits. The 
statistical errors and the covariance matrix were deter- 
mined with the jackknife method. The fit ranges were 
chosen such that they extend over those values of t where 
we see a credible effective mass plateau. Changing the 
upper limit of the fit interval by ±1 does not affect the fit 
results (the variation is considerably less than one s.d.). 
Also an increase of the lower limit of the fit range changes 
the fit results by considerably less than one s.d. When 
decreasing the lower end of the fit interval, the result for 
the mass typically goes up by one s.d. This effect is obvi- 
ous from the effective mass plots where one sees that for 
smaller t one runs into the contributions of higher excited 
states. Our fits all have \ 2 /d.o.f. values smaller than 1 
and the % 2 /c?.o./. also stays below 1 when considering 
the changes in the fit range. Thus, the x 2 /d.o.f. is not a 
very stringent criterion for determining the fit range and 
our choice based on the effective mass plateaus is rather 
conservative. 

In our fits we carefully analyzed the possibility of 
normalizing our matrix as suggested by the generalized 
eigenvalue problem Eq. @. In particular we tried the 
choices to — 0,1,2 and found that the fits for the masses 
are unchanged within error bars. For the ground state 
the change was less than 0.3 s.d. and for the excited state 
less than 1 s.d. For to = 2 we find a noticeable decrease 



in the quality of the effective mass plateaus, which comes 
from the fact the at <o = 2 the matrix C(io) -1 which is 
used for the normalization already becomes affected by 
statistical noise. We cannot confirm any suppression of 
the effects from higher excited states and our results from 
the generalized eigenvalue problem are indistinguishable 
from the data obtained from the regular eigenvalue prob- 
lem. In the further discussion we use the latter. 

We experimented also with another possibility to fit 
our correlators. Based on the Hilbert space decomposi- 
tion we used an ansatz of the form 



Cij(t) = 



R 

E 

n=l 



A (n) A ( r ] e- 



L M„. 



(17) 



to fit an R x R correlation matrix. The coefficients 
are chosen real, since the correlation matrix is real and 
symmetric. We find that the ansatz 1|17[) can be applied 
successfully only for 2x2 matrices. For larger correlation 
matrices the data of our ensembles used for testing the 
method are not accurate enough for the multi-parameter 
fit i|17|) . The x 2 /d.o.f. becomes large and we sometimes 
encountered convergence problems in the numerical mini- 
mization of the x 2 functional. For the 2x2 case, however, 
the results from fitting the whole correlation matrix with 
(|17fl agree reasonably well with the fits of the eigenval- 
ues. For the ground state mass the two results differ by 
less than one s.d. For the first excited state the differ- 
ence is below one s.d. for am — 0.04 rising to 1.5 s.d. at 
am = 0.20. 

In Fig. 0] we show the fitted nucleon masses in lattice 
units as a function of the quark mass. Open triangles rep- 
resent our data and the statistical error for the ground 
state and the filled triangle gives the experimental mass 
for the nucleon iV(940) (converted to lattice units using 
the Sommer parameter). The open circles are the data 
for the first excited state and the corresponding filled cir- 
cle represents the mass of the Roper iV(1440). The open 
squares are the masses for the lowest negative parity state 
calculated on the same ensemble of configurations in Q . 
Note that the lower two data sets are for ground states of 
positive and negative parity and thus have smaller sta- 
tistical errors. Although we do not attempt a chiral ex- 
trapolation of our data (this is deferred to a future large 
scale study of the nucleon system), our numbers seem 
to approach the experimental data reasonably well. For 
the quark mass range in our study the first excited, posi- 
tive parity masses still remain above the negative parity, 
ground state data, but a trend towards a level crossing 
is plausible. We remark that the smallest quark mass 
where we can fit the excited nucleon, am = 0.04, cor- 
responds to a pion mass of 450 MeV (see 0). Thus 
our preliminary data contradict a level crossing at 600 
MeV as claimed in Clearly further studies in larger 
volumes are important; such studies are in progress. 
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IV. EXAMPLE B: EXCITED RHO MESON 

As another test of our approach we discuss the rho- 
meson p(770) and its radial excitation p(1450). Tradi- 
tionally the /o(770) has been studied with the vector cur- 
rent interpolator, 



Eigenvalues 



Effective masses 



u(x) 7j d(x) . 



(18) 



The same meson in the chiral symmetry broken regime 
should also be seen with the time component of the tensor 
interpolator, 



u(x) an d(x) 



(19) 



The difference between the two is that they transform 
under different representations with respect to SU(2)l X 
SU(2)n and U(1)a Wh - Again we use wide and nar- 
row quark sources for both interpolators. Thus, for both 
interpolators (JTSJ, l|19fl we can build the 4 operators as 
listed in Eq. I|10|) . However, the two combinations (n,w) 
and (w, n) give identical correlators and one of them can 
be omitted. Thus, we evaluate two 3x3 correlation ma- 
trices, one for the vector interpolator and one for the 
tensor interpolator (the different spatial components are 
averaged). For both these interpolators we check whether 
the ground and radially excited states couple. When we 
diagonalize the 3x3 matrix with either the interpolator 
(|18|l or the interpolator l|19fl we see a pronounced ex- 
ponential decay only for the two larger (in magnitude) 
eigenvalues, X^\t) and \^ 2 \t). The smallest eigenvalue 
A^ 3 ) (t) does not show a clear effective mass plateau and 
becomes negative at small quark masses for large t. This 
is a clear indication that this eigenvalue couples to an 
unphysical quenched ghost state |23tl24|. 

In Fig. we show the two larger eigenvalues and the 
corresponding effective mass plots obtained with the ten- 
sor current (|19fl for three different quark masses (am — 
0.05, 0.10, 0.20, from top to bottom). We observe a clean 
exponential decay for both eigenvalues and correspond- 
ing plateaus in the effective mass plots. We remark, that 
results for the mass obtained with the vector current i|18|) 
agree within the error bars. Contrary to the baryon case, 
we observe only one radial excitation. Indeed, there is 
only a single one-node excitation for a meson with the 
given quantum numbers. 

Like for the baryons, we apply a standard correlated 
single exponential fit to extract the meson masses. For 
the ground state the time interval is chosen to be (3,8), 
while for the excited state we use (2,5). The final results 
for the masses as a function of the quark mass are shown 
in Fig. [SI We find that the ground state approaches its ex- 
perimental value reasonably well (the experimental data 
were converted to lattice units with the Sommer param- 
eter scale). 

The excited state masses are considerably above their 
experimental value. There is, however, a plausible reason 
for this behavior. The sizes of hadrons which are not, 
or only weakly, affected by spontaneous chiral symmetry 
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FIG. 5: Left hand side panels: The largest eigenvalue (circles) 
and the next to the largest one (squares). Right hand side: 
The corresponding effective mass plots. From top to bottom 
the quark masses are am — 0.05, 0.10, 0.20. 



breaking can be estimated from the known string tension, 
which is approximately 1 GeV/fm. Hence the size of the 
p-meson is expected to be below 1 fm, while the size 
of p(1450) should be approximately 1.5 fm. The size of 
our lattice is 1.8 fm, which is clearly not enough for a 
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precise measurement of the p(1450) mass. The finite size 
effect cannot be neglected for the excited state and shifts 
the measured mass up as compared to the experimental 
value. A study of the rho system on larger lattices is in 
progress. 



V. ANALYZING THE OPERATOR CONTENT 
OF THE PHYSICAL STATES 

We have based our new approach on the working hy- 
pothesis that, when using the variational method, the 
excited states have a better overlap with a basis of oper- 
ators that allow for a node in the spatial wave function. A 
crucial test of this assumption is to check whether indeed 
the ground state is built from a nodeless combination of 
our sources and the excited states do show nodes. 

This question can be addressed by analyzing the eigen- 
vectors of the correlation matrix JHJ). Let us denote by 
e w the fc-th eigenvector of the correlation matrix, and 
its eigenvalue is A^. Then we can define optimal oper- 
ators Ok by 



O k 



E 



c {k) 0- 



(20) 



with the mixing coefficients Cj determined from the en- 
tries of the eigenvector e^ k > via 



„(*) 



(21) 



where the asterisk denotes complex conjugation. The 
correlation matrix of the optimal operators Ok is diago- 
nal, i.e. 



(O fc (t)Oz(0)t) = fc,A«(t) . 



(22) 



Thus, the operator 0\ has optimal overlap with the 
ground state, the operator O2 has optimal overlap with 
the first excited state etc. By analyzing the coefficients 

(k) 

Cj we can thus learn about the structure of the state 

seen in the eigenvalue . In particular we can address 
the question whether there are nodes in the wave func- 
tions. As we have discussed in Section 2, a node occurs 
if there is a relative minus sign between the wide and the 
narrow quark source. For the nucleon system, where we 
work with the set of operators listed in i|15|) . our work- 
ing hypothesis implies that all coefficients have the same 
sign for the ground state, while for the excited states we 
expect relative signs. 

We determine the eigenvectors of C(t) at each 
value of t. Since C(t) is real and symmetric, the eigenvec- 
tors can be chosen real. In Fig.[7|we display these eigen- 
vectors for the nucleon analysis by plotting the (real) 

coefficients = . In the top plot we show the num- 
bers for the ground state (k = 1), in the middle plot the 
numbers for the first excited state (k = 2) and in the bot- 
tom plot we display the second excited state (k = 3). The 
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FIG. 7: Entries of the eigenvectors as a function of t. The 
data are for quark mass am = 0.10. The entries are labeled 
by the source combinations they correspond to. 



different symbols represent the mixing coefficients for the 
different basis operators and we use circles for (n,n,n), 
triangles for (n, n, w), upside down triangles for (n, w, n) 
and squares for (w, w, n), all referring to the interpolator 
Xi- The data we show are for quark mass am = 0.10. 
We remark that our eigenvectors are normalized to one, 
implying that 



J»2 



+ c 



(fe)2 



= 1, 



(23) 



for the mixing coefficients. Furthermore, since C(t) is 
real symmetric, the eigenvectors are also orthogonal to 
each other for all t. 

It is obvious from the plot that also the mixing coeffi- 
cients show plateaus. For the ground state these plateaus 
start at t = 0, while for the excited states they typi- 
cally start at t = 4, where we also observed the onset 
of plateau-like behavior in the effective masses. Starting 
at t = 8 the plateaus for the coefficients vanish for the 
excited states. At this t the eigenvalues are already very 
small (they decay exponentially) and different eigenvec- 
tors start to mix. However, for all quark masses we ob- 
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FIG. 8: Coefficients for the optimal operator combinations as 
a function of the quark mass. 



serve long enough plateaus to clearly identify the mixing 
coefficients of the optimal operators. 

In Fig.|S|we show the mixing coefficients as a function 
of the quark mass. In particular we choose the values 
of the coefficients for time slice t = 5. We use the same 
symbols for the different basis operators as in Fig. and 
again show from top to bottom the ground state, the first 
and the second excited state. 

The top plot shows that for the ground state (the nu- 
cleon) all 4 mixing coefficients have the same sign and the 
wave function of the ground state does not have a node. 
For the two excited states the situation is clearly differ- 
ent. We find positive and negative coefficients, indicating 
that the wave function has a node. We stress that a lin- 
ear combination of the profiles shown in Fig. ^ using the 
coefficients of Fig.[S]should not be literally interpreted as 
the true nucleon wave function. Firstly, our narrow and 
wide sources are not orthonormalized, and secondly they 
provide a much too small basis for mapping the details 
of a complicated three body wave function. 

We performed the same analysis also for the eigenvec- 
tors of the correlation matrices we used for the rho me- 



son. Again we confirm that the ground state has equal 
sign coefficients and thus is nodeless, while the excited 
state has relative signs between its coefficients indicat- 
ing a node. Thus for both systems where we tested our 
method we could confirm that in the variational method 
combinations without node give rise to the ground state 
signal while the signal for excited states comes from com- 
binations with a node. 



VI. SUMMARY 

We have presented a new approach to improving the 
spatial structure of operators used for hadron spec- 
troscopy on the lattice. The central idea is to combine 
Jacobi smeared quark sources of different width in the 
variational approach. This combination allows for nodes 
in the spatial wave function and improves the signal from 
radially excited states. We have demonstrated the power 
of this approach by applying the method to a study of 
the lowest radial excitation of the nucleon, the Roper 
resonance, and to the radial excitation of the rho me- 
son, p(1450). In both cases we clearly identify credible 
plateaus in the effective mass plots and are able to trace 
the signal of these states from the heavy quark region 
towards the chiral limit. It is reassuring to see clear 
effective mass plateaus also for the excited states, and 
the corresponding mass needs not be extracted from a 
multi-exponential fit. A simple single exponential fit is 
sufficient for extracting the mass from the decay of the 
second eigenvalue. When plotting the fitted masses as a 
function of the quark mass, the agreement with the ex- 
perimental situation is reasonable. A large scale study 
with larger lattices and a systematic chiral extrapolation 
is in progress. We also note that it is possible to ex- 
tend our method to the next radial excitations, with two 
nodes. For that one needs to use at least three different 
types of quark sources. We plan to study this extension 
of the method in the near future. 

There are two physical implications of our preliminary 
observations: 

(i) The fact that our numerical data for the first excited 
state can be plausibly interpreted as the Roper signal has 
an important consequence for the nature of this state: 
It implies that the Roper's leading Fock component is 
a 3 quark state. A trend towards a level switching of 
the lowest positive and negative parity excited states is 
visible, which is consistent with the physical picture of 
Ref. [Hi. 

(ii) In our data the quark mass dependence of the ex- 
cited rho meson state p(1450) is only weak. This im- 
plies that p(1450) is only weakly affected by spontaneous 
breaking of chiral symmetry, in contrast to the lowest 
baryon states. This indicates that the physics of the 
excited meson states and of the lowest baryons is very 
different and that perhaps we observe a transition to the 
restoration of chiral symmetry in highly excited hadrons 
ii. 
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These are important conclusions and we expect that 
the application of our source techniques in a large scale 
study can lead to a deeper understanding of these issues 
from a lattice perspective. 
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